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A  hybrid  system  mainly  composed  of  a  DCFC  (Direct  Carbon  Fuel  Cell),  a  TEG  (Thermoelectric  Generator) 
and  a  regenerator  is  put  forward,  where  the  DCFC  electrochemically  converts  the  chemical  energy  in  the 
solid  carbon  into  electricity  and  waste  heat,  and  the  TEG  further  converts  the  waste  heat  into  electricity 
for  additional  power  generation.  The  main  irreversibilities  in  each  element  of  the  hybrid  system  are  char¬ 
acterized,  and  the  heat  losses  between  the  DCFC  and  the  environment  are  also  considered.  Numerical 
expressions  for  the  power  output  and  efficiency  of  the  hybrid  system  are  respectively  derived,  from 
which  the  general  performance  characteristics  are  revealed.  The  fundamental  relationship  between  the 
operating  current  density  of  the  DCFC  and  the  dimensionless  electrical  current  of  the  TEG  is  obtained, 
and  thus  the  region  of  the  operating  current  density  of  the  DCFC  that  the  TEG  exerts  its  function  is  deter¬ 
mined.  By  using  such  a  system,  the  equivalent  maximum  power  density  of  the  hybrid  system  allows  to  be 
50%  larger  than  that  of  the  sole  DCFC  system.  The  effects  of  the  operating  current  density,  operating  tem¬ 
perature,  heat  conductivity,  and  some  integrated  parameters  on  the  performance  of  the  hybrid  system 
are  discussed. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  worldwide  growing  demand  for  electrical  energy  cannot  be 
compensated  by  renewable  energies  in  the  medium  term  [1],  Coal 
is  the  most  abundant  and  economic  fossil  resource  on  the  earth 
and  is  forecast  to  account  for  28%  of  the  word’s  increased  primary 
energy  consumption  by  2030  [2],  Therefore,  the  development  and 
deployment  of  clean  coal  technologies  are  crucial  to  promote  sus¬ 
tainable  development  [3],  A  promising  development  in  this  context 
could  be  direct  carbon  fuel  cell  (DCFC)  that  has  the  potential  to 
directly  convert  the  chemical  energy  in  solid  carbon  into  electricity 
without  the  need  for  gasification  or  the  moving  machinery  associ¬ 
ated  with  conventional  electric  generators.  DCFCs  offer  compara¬ 
tive  high  efficiency  of  conversion  and  can  yield  2-3  times  the 
amount  of  electrical  energy  as  coal-fired  power  plants  for  a  given 
amount  of  carbon,  resulting  in  a  50%  reduction  in  carbon  dioxide 
emissions  per  unit  of  electricity  generation  [4],  DCFCs  can  use  coal 
coke,  petroleum  coke,  and  gas  carbon  as  well  as  biomass  carbon  as 
fuel,  and  the  exhaust  gas  is  an  almost  C02  stream  [5,6],  The  nearly 
pure  C02  stream  can  be  easily  collected  for  downstream  disposal 
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using  methods  such  as  carbon  sequestration,  and  thereby,  the 
energy  and  cost  penalties  to  capture  the  C02  will  also  be  signifi¬ 
cantly  less  than  other  technologies  [7], 

However,  one  of  the  major  drawbacks  of  DCFCs  that  prohibits 
its  commercialization  stage  is  its  comparative  low  power  density 
[7],  The  maximum  power  density  obtained  from  DCFCs  has  under¬ 
gone  many  advances  in  recent  years  not  only  due  to  the  progresses 
in  the  materials  fabrication  but  also  due  to  the  optimized  operating 
conditions  and  DCFC  system  design  [8-14],  As  an  alternative 
approach,  the  equivalent  maximum  power  density  of  DCFC  can 
be  also  effectively  elevated  by  employing  cogeneration  systems 
[15-17],  Steinberg  [15]  developed  a  high  efficiency  energy  cycle 
by  integrating  a  molten  carbonate  DCFC  with  a  hydrogen  plasma 
black  reactor  (HPBR)  to  convert  the  fossil  and  biomass  fuels  to  elec¬ 
tric  power  and  transportation  fuels,  it  was  found  that  the  overall 
thermal  efficiency  for  electricity  and  hydrogen  or  transportation 
fuel  was  estimated  to  vary  from  70%  to  90%.  Liu  et  al.  [16]  built 
and  simulated  a  highly  efficient  energy  system  mainly  composed 
of  a  methane  catalytic  decomposition  reactor  (MCDR),  a  DCFC,  an 
internal  reforming  solid  oxide  fuel  cell  (1RSOFC)  and  two  gas  tur¬ 
bines.  They  found  that  the  DCFC  had  better  high-load  performance, 
while  the  IRSOFC  performed  better  at  relatively  lower  loads, 
and  80%  of  the  C02  emission  was  reduced  compared  to  a  system 


utilizing  a  single  SOFC.  Their  further  studies  showed  that  the  elec¬ 
trical  exergy  efficiency  of  68.24%  for  the  system  was  achieved  and 
the  combined  power-heat  exergy  efficiency  could  be  over  80%  if 
the  waste  heat  in  the  exhaust  was  recovered  [17], 

TEG  (Thermoelectric  Generator)  can  directly  convert  thermal 
energy  into  electricity  via  the  Seebeck  effect;  that  is,  when  a  tem¬ 
perature  difference  is  created  across  the  thermoelectric  device  and 
an  external  load  is  properly  connected,  a  direct-current  voltage 
develops  in  the  thermoelectric  elements  and  electrical  current 
flows  through  the  load  [18,19],  Compared  with  conventional  elec¬ 
tric  power  generators,  TEG  offers  many  advantages  such  as  being 
highly  reliable,  having  no  moving  parts,  and  being  environmentally 
friendly  [20],  There  is  a  growing  interest  for  waste  heat  recovery 
TEG  using  various  heat  sources  such  as  combustion  of  solid  waste, 
geothermal  energy,  power  plants,  and  other  industrial  heat-gener¬ 
ating  processes  [21  -23],  When  a  DCFC  is  put  into  operation,  a  great 
part  of  chemical  energy  in  the  fuel  will  be  finally  released  as  waste 
heat  due  to  the  Joule  effect  of  the  equivalent  internal  resistance. 
Obviously,  the  waste  heat  generated  in  the  DCFC  can  be  readily 
converted  into  power  generation  through  a  bottoming  TEG  and 
thus  the  performance  of  the  DCFC  can  be  effectively  enhanced. 

In  this  work,  a  hybrid  system  consisting  of  a  DCFC,  a  TEG  and  a 
regenerator  is  proposed  to  recover  the  waste  heat  generated  in  the 
DCFC.  The  performance  parameters  of  the  hybrid  system  are  ana¬ 
lytically  derived  by  considering  the  main  irreversible  losses  in  each 
element  of  the  system.  The  effectiveness  of  the  hybrid  system  will 
be  shown  by  some  numerical  calculation  examples,  and  the  effects 
of  some  operating  conditions  and  designing  parameters  on  the  per¬ 
formance  of  the  hybrid  system  will  be  revealed. 


2.  A  DCFC-TEG  hybrid  system 

The  hybrid  system  presented  here  consists  of  a  DCFC,  a  TEG  and 
a  regenerator,  as  schematically  shown  in  Fig.  1,  where  the  TEG  is 
operated  between  the  heat  source  (i.e.,  the  DCFC)  at  temperature 
T  and  the  heat  sink  (i.e.,  the  environment)  at  temperature  T0,  qi 
and  q2  are  respectively  the  rate  of  heat  input  from  the  DCFC  to 
the  thermoelectric  generator  and  the  rate  of  heat  rejection  from 
the  thermoelectric  generator  to  the  environment,  qL  is  the  rate  of 
heat  loss  between  the  DCFC  and  the  environment,  Pdcfc  and  Ptec 
are  respectively  the  electrical  power  outputs  of  the  DCFC  and 
TEG.  In  this  hybrid  system,  the  DCFC  not  only  works  as  an  electro¬ 
chemical  power  generator  but  also  acts  as  the  high-temperature 
heat  reservoir  of  the  TEG  which  generates  some  additional  electric¬ 
ity,  and  the  regenerator  is  applied  to  preheat  the  incoming  fuel  and 
oxidant  with  the  comparative  high-temperature  exhaust  product 
from  the  DCFC.  Employing  such  a  hybrid  system,  the  heat  produced 


in  the  DCFC  can  be  effectively  harvested,  and  consequently,  the 
performance  of  the  DCFC  can  be  improved. 

The  whole  hybrid  system  is  formulated  based  on  the  following 
assumptions  [24-28]: 


(1 )  Both  the  DCFC  and  the  TEG  are  operated  under  steady-state 
conditions. 

(2)  Geometric  configuration  of  the  TEG  is  in  the  optimum  form. 

(3)  Operating  temperature  and  pressure  are  uniform  and  con¬ 
stant  in  the  DCFC. 

(4)  Complete  chemical  reactions  are  considered  and  no  reactant 
is  left  after  the  reactions  [27], 

(5)  All  gases  involved  are  assumed  to  be  compressible  ideal 
gases. 

(6)  Carbon  fuel  is  regarded  as  a  rigid  sphere  and  packed  with  a 
simple  hexagonal  pattern. 

(7)  Carbon  fuel  is  assumed  to  be  oxidized  by  completely  electro¬ 
chemical  reactions. 

(8)  Heat  leakage  through  the  TEG  surroundings  is  neglected. 

Based  on  the  above  simplifications  and  assumptions,  each  com¬ 
ponent  in  the  hybrid  system  will  be  analyzed,  and  then  the  perfor¬ 
mance  characteristics  of  the  hybrid  system  will  be  synthetically 
investigated. 


2.J.  The  DCFC 


The  DCFC  directly  converts  the  chemical  energy  of  the  incoming 
graphite  into  electricity  and  heat  without  any  gasification,  it  is 
mainly  composed  of  a  packed  bed  anode  and  a  lithium  doped 
NiO  cathode  with  U2CO3/K2CO3  eutectic  melt  as  an  electrolyte 
sandwiched  between  the  two  electrodes.  The  overall  electrochem¬ 
ical  reaction  proceeded  in  a  DCFC  can  be  summarized  as  C  +  02  -» 
C02  +  Electricity  +  Heat,  electrons  are  drawn  from  the  anode  to  the 
cathode  through  an  external  circuit  and  a  direct  electric  current  is 
produced.  As  previously  described  in  Ref.  [25],  the  output  cell 
voltage  of  a  DCFC,  V,  is  always  lower  than  the  OCV  (open  circuit 
voltage),  E,  because  there  are  some  voltage  losses  resulting  from 
the  electrode  kinetics,  electronic/ionic  resistances  of  the  cell  com¬ 
ponents,  and  reactants/product  transportation,  which  correspond 
to  the  activation  overpotential  ( Vact ),  ohmic  overpotential  {Vohm), 
and  concentration  overpotential  (Vcon),  respectively.  Taking  these 
three  losses  into  account,  one  may  obtain  the  output  cell  voltage, 
power  output,  and  efficiency  of  a  DCFC  at  certain  operating  current 
density 


V  =  E 
Pdcfc  =  VI  =  VJA, 


Vaa 


-Vcon 


-vohn 


(1) 

(2) 


neFV 
Ah  ’ 


(3) 


where  E  =  £„  +  §  In  [po2mt(Pco2cat)2Pci2,an]  is  the  OCV,  E0  =  -Ag° 
(T)/(neF),  £0  and  Ag°(r)  are,  respectively,  the  reversible  potential 
and  the  Gibbs  free  energy  change  at  the  standard  pressure  (1  atm) 
and  temperature  T,  ne  is  the  number  of  electrons  involved  per  reac¬ 
tion,  R  is  the  universal  gas  constant,  pCh  and  pC02  are,  respectively, 
the  partial  pressures  of  02  and  CO2  at  the  electrode  surfaces,  and 
the  subscripts  “an"  and  “cat”  represent  anode  and  cathode,  respec¬ 


tively;  Vc 


t,n,  =  i  In  {//  (2Jo,m)  +  J\J/(2J0,m)}2  4 


s  the  activation 


1,  /  and  ]  are  respectively  the  el 
ating  current  density  of  the  DCFC,  the  si 
or  "cat”,  J0,an  =  Kb  exp  (-£b/T)  and  J0|COt  -  fBcat(pC02,atY'  (Po^atY2  are 
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respectively  the  anode  exchange  current  density  and  the  cathode 
exchange  current  density  [29,30];  Vcon  =  §  ln[/lim/(/lim  - J)]  is  the 
concentration  overpotential,  Jlim  =  neFKc o2  CbCOi  is  the  limiting 
current  density,  KCo2  is  the  mass  transportation  coefficient  of  C02, 
and  Cbc 02  is  the  concentration  of  C02  in  the  bulk  [24,31];  Vohm  = 
X  IcjRcj  +  X  h.iHej  +  14c  is  the  ohmic  overpotential,  /c,  and  lei  are 
respectively  the  electric  currents  flowing  through  the  carbon  phase 
and  electrolyte  phase  at  the  ith  slab  in  the  anode,  RCj  and  Rei 
are  respectively  the  resistances  in  the  carbon  phase  and  electrolyte 
phase  at  the  ith  slab  in  the  anode,  Vec  is  the  total  ohmic  overpotenial 
in  the  electrolyte  and  the  cathode  [25];  -AH  =  -  0/A  is  the  total 
chemical  energy  released  per  unit  time,  (-Ah)  is  the  molar  enthalpy 
change  of  the  electrochemical  reaction  [25],  A  is  the  polar  plate  area 
of  the  DCFC. 


2.2.  The  regenerator 

The  regenerator  in  the  hybrid  system  functions  as  a  counter¬ 
flow  heat  exchanger,  which  heats  the  incoming  reactants  from 
the  ambient  temperature  to  the  operating  temperature  of  DCFC 
by  absorbing  the  comparative  high-temperature  heat  contained 
in  the  product.  Owing  to  the  existence  of  thermal  resistance,  some 
regenerative  losses  are  inevitable  and  some  additional  heat  from 
the  DCFC  should  be  transferred  to  compensate  the  regenerative 
losses  immediately.  Assuming  the  rate  of  the  regenerative  losses 
is  proportional  to  the  temperature  difference  between  the  DCFC 
and  the  environment,  the  regenerative  losses  is  given  by  [32,33] 

qre  =  areAre(l-£)(T-To),  (4) 

where  are  and  Are  are  the  heat-transfer  coefficient  and  heat-transfer 
area  of  the  regenerator,  respectively. 


2.3.  The  TEC 


As  show  in  Fig.  2,  the  TEG  in  the  hybrid  system  is  closely 
attached  to  the  DCFC,  it  is  composed  of  many  n-  and  p-type  semi¬ 
conductor  legs  that  are  connected  electrically  in  series  by  metal 
strips  and  thermally  in  parallel,  the  thermoelectric  generating  ele¬ 
ments  are  assumed  to  be  insulated,  both  electrically  and  thermally, 
from  their  surroundings,  except  at  the  junction-reservoir  contacts. 
Ig  is  the  current  assumed  to  flow  along  the  arms  of  the  generator,  T, 
and  T2  are  the  temperatures  of  the  hot  junction  and  cold  junction 
respectively,  for  simplification,  Tt  is  assumed  to  be  equal  to  the 
operating  temperature  of  the  DCFC,  T,  and  the  temperature  of  the 
cold  junction,  T2,  is  assumed  to  be  equal  to  the  temperature  of 
the  environment,  T0.  When  the  TEG  operates  stably,  the  boundary 
conditions  are  determined  by  Ti(0)  =  T2(0)=-  ■  T,<0)  =  ■  ■  ■  =  Tm( 0)  =  T 

and  T1(L)  =  T2(L)=  -  T,{L) - Tm(L)  =  T0,  where  L  is  the  arm 

length  of  the  TEG  and  m  is  the  number  of  TEG  couples.  Owing  to 
the  fact  that  there  exists  a  finite  temperature  difference  between 
the  hot  and  the  cold  junctions,  a  heat  leak,  K(T-T0),  from  the  hot 
junction  to  the  cold  junction  via  the  semiconductor  elements 
should  be  considered  [34],  Thus,  the  heat  conductions  of  the  TEG 


between  the  DCFC  and  the  environment  can  be,  respectively, 
expressed  as  follows  [35,36] 


qx  =  aIgT  -  0.5 12gR  +  K(T  -  T0), 

(5) 

and 

q2  =  (XlgT0  +  0.5/gR  +  K(T  -  T0), 

(6) 

where  a  =  (ap  -  a „)m  is  the  total  Seebeck  coefficient,  the  subscripts 
“p"  and  “n”  stand  for  p-  and  n-type  semiconductor  legs;  R  =  [pplpt 
SP  +  pJnlS„)m  is  the  total  electrical  resistance  of  the  semiconductor 
couple,  p  is  the  electrical  resistivity,  l  is  the  length  of  the  semicon¬ 
ductor  arms,  S  is  the  cross-sectional  areas  of  semiconductor  arms; 
IgR  represents  the  Joule  heat  resulted  from  the  electrical  resistance 
in  the  TEG;  K  =  {Kjfipfip  +  KnSn//„)m  is  the  total  thermal  conductance 
of  semiconductor  couple,  k  is  the  thermally  conductivity  of  the 
semiconductor  materials. 

Based  on  Eqs.  (5)  and  (6),  the  power  output  and  efficiency  of  the 


TEG  can  be,  respectively,  expressed  as 

Ptec  =  q\~q2  =  «(T  -  T0)j  -  kj2 /Z, 

(7) 

and 

7(7  r°)J  , 
q,  Z(T  -  T0)  +  ZTj  -  J2/2 

(8) 

where  j  =  odg/K  and  Z  =  a2/(KR)  are,  respectively,  the  dimensionless 
electrical  current  and  the  figure  of  merit  of  the  semiconductor 
materials  [37],  When  some  heat  are  transferred  from  the  DCFC 
DCFC  to  the  TEG,  both  the  power  output  and  the  efficiency  are  lar¬ 
ger  than  zero,  i.e.,  Ptec  >  0  and  rime,  >  0,  and  thus  one  may  easily 
determine  the  region  of  the  dimensionless  electrical  current  that 
the  TEG  takes  effect  as  follows 

0<j<Z(T-T0),  (9) 

where  Z(T  -  T0 )  is  the  maximum  dimensionless  electrical  current 
allowable  for  the  TEG  to  generate  electricity. 

2.4.  The  efficiency  and  power  output  of  the  hybrid  system 


As  clearly  illustrated  in  Fig.  1,  the  total  heat  released  in  the 
DCFC  is  divided  into  three  parts.  Besides  the  part  used  to  replenish 
the  regenerative  losses  as  described  by  Eq.  (4),  the  part  leaked  into 
the  environment  via  convective  and/or  conductive  heat  transfer 
and  the  part  transferred  to  the  TEG  for  electricity  generation  can 
be  respectively  expressed  as 

qL  =  aLAL(T-T0),  (10) 


neFcfiT  -  T0)  neFc2(T-T 0)1 

-Ah  J  ■ 


where  is  the  convective  and/or  conductive  heat-leak  coefficient, 
Al  is  the  heat-transfer  area,  c,  =  ar<Are(l  -  e)/A  and  c2  =  a,A,/A  are 
two  constants  which  are  independent  of  temperature. 

By  considering  the  heat  leak  between  the  hot  and  cold  junc¬ 
tions,  the  TEG  in  the  hybrid  system  begins  to  exert  its  function  only 
when  the  following  condition  is  satisfied: 

-AH  -  PDcfc  >  qre  +  Ql  +  K(T  -  T0).  (12) 

Substituting  Eqs.  (2),  (4)  and  (10)  into  Eqs.  (11),  (12)  can  be  fur¬ 
ther  explicitly  rewritten  as 

> -  [=gfSsd l(Ci - cMT  ~  T") + *(T  -  ™  (,3) 


Fig.  2.  The  schematic  diagram  of ; 


jple  TEG. 
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Parameters  used  in  the  modeling  [24,25,37,29,38], 


Parameter  Value 

Operating  pressure,  P  (Pa) 

1.01325  x  105 

Ideal  gas  constant,  R  (J/mol  K) 

8.314 

Faraday  constant,  F  (C/mol) 

96,485 

Height  of  packed  bed  anode,  H  (m) 

1.0  x  10“3 

Number  of  electrons  involved  per  reaction,  ne 

4 

Polar  plate  area  of  a  DCFC,  A  (m2) 

0.04 

Diameter  of  spherical  graphite  particle,  Dc  (m) 

1.0  x  10~5 

[24,25] 

Concentration  independent  cathode  exchange  current 

5.0  x  102  [25] 

density,  jo  cat  (A/m2) 

Pre-exponential  factor  of  the  backward  reaction,  KB  (A/m2) 

5.8  x  109  [29] 

Temperature  activation  of  the  backward  reaction,  EB  ( K  1 ) 

22,175  [29] 

Mass  transport  coefficient  of  C02,  Kc o2  (m/s) 

3.5  x  10~2  [25] 

Constant,  r, 

-1.250  [24,25] 

Constant,  r2 

0.375  [25] 

Bulk  gas  compositions  at  the  cathode  side 

67%C02/33%  02 
[25] 

Temperature  of  environment,  T0  (K) 

298 

Constant  c,  (W  m  2  K  ’) 

0.1 

Constant  c2  (W  m  2  K  ’) 

0.1 

Figure  of  merit  of  the  semiconductor  materials,  ZT0 

1.0  [37] 

Heat  conductivity,  K  (W  K  1  m  ’) 

1.5  x  10  2  [38] 

3.  Results  and  discussion 

As  shown  by  Eqs.  (16)-(19),  the  performance  of  the  TEG-based 
DCFC  hybrid  system  not  only  depends  on  the  parameters  related  to 
the  DCFC  such  as  operating  temperature,  operating  current  density, 
electrolyte  thickness,  but  also  depends  on  the  parameters  related 
to  the  TEG  such  as  electrical  resistance  and  thermal  conductance, 
figure  of  merit  of  the  thermoelectric  materials,  geometric  parame¬ 
ters  and  so  on.  In  the  following,  numerical  calculations  are  carried 
out  by  using  the  commercial  software  MATLAB®  based  on  the 
mathematical  model  developed  in  Section  2  and  the  parameters 
summarized  in  Table  1,  and  these  parameters  are  taken  as  default 
values  unless  specifically  mentioned. 

The  power  densities  and  efficiencies  of  the  DCFC,  TEG  and 
hybrid  system  varying  with  the  operating  current  density  of  the 
DCFC  are  clearly  shown  in  Fig.  3,  where  P*  =  P/A  is  the  power  den¬ 
sity,  is  the  maximum  power  density,  and  JP  is  the  operating 
current  density  corresponding  to  P*max.  Fig.  3(a)  shows  that  the 
state  of  the  hybrid  system  at  the  maximum  power  output  is  differ¬ 
ent  from  that  of  the  DCFC  or  TEG  at  the  maximum  power  output.  It 
is  also  observed  that  the  power  density  of  the  hybrid  system  first 
increases  and  then  decreases  as  the  operating  current  density  is 
increased,  while  the  efficiency  first  rapidly  decreases  then  slightly 


where  Jc  is  the  critical  operating  current  density  of  the  DCFC  from 
which  the  TEG  in  the  hybrid  system  starts  to  work. 

When  J  >  Jc  and  0  <j  <Z(T-  T0),  the  fundamental  relationship 
between  the  operating  current  density  of  the  DCFC,  J,  and  the 
dimensionless  current  of  the  TEG,  j,  is  given  by  the  following 
equation: 


f  ~  2ZTj  -  2 Z(T  -ro)  +  2,f  (1  -  r,DCFC)J  -  (ct  +  C2)(T  -  T„)] 

-0.  (14) 

Combining  Eqs.  (9)  and  (14),  one  may  easily  determine  the 
maximum  operating  current  density  of  the  DCFC,  JM,  from  which 
the  TEG  in  the  hybrid  system  stops  working.  Thus,  the  TEG  in  the 
hybrid  system  cannot  take  effects  in  the  entire  operating  current 
density  region  of  DCFC  but  only  works  in  the  following  region 

Jc<J<Jm ■  (15) 

When JC<J< Jmi  based  on  Eqs.  (2),  (3),  (7)  and  (8),  the  numerical 
expressions  of  the  power  output  and  efficiency  of  the  hybrid  sys¬ 
tem  can  be,  respectively,  expressed  as 

P  =  Pdcfc  +  Pteg  =  VJA  +  I<T0[(T/T0  -  1  )j  -  f/(ZT0)} ,  (16) 

T1  Pdcfc  +  Pteg  „  ,  neFKT0[(T/T0  -  l)j  -f/(ZT0)} 

n  w  ^/aaS  ' 

When  J  ^Jc  or  J  >  JM,  the  power  output  and  efficiency  of  the 
hybrid  system  equal  to  that  of  the  sole  DCFC,  i.e., 

P  =  Pdcfc,  (18) 

and 


V  —  tIdcfc-  (19) 

It  is  worthwhile  to  mention  that  J  and  j  in  Eqs.  (16)  and  (17)  are 
not  independent  of  each  other,  and  the  relationship  between  J  and  j 
is  given  by  Eq.  (14). 


Fig.  3.  The  curves  of  (a)  power  densities  and  (b)  efficiencies  of  the  DCFC,  TEG,  and 
hybrid  system  varying  with  the  operating  current  density  of  the  DCFC,  where  P*  =  P/ 
A  is  the  power  density  of  the  hybrid  system,  P*DCFC  and  P^c  are  the  power  densities  of 
the  DCFC  and  TEG ,JP  is  the  current  density  corresponding  to  the  maximum  power 
density  of  the  hybrid  system  PJ^. 
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increases  and  finally  decreases  as  the  operating  current  density 
increases.  The  intermediate  slightly  augment  in  the  efficiency 
occurs  in  the  region  of  Jc  < J  <  Jm  because  the  increase  in  the  effi¬ 
ciency  contributed  by  the  TEG  is  larger  than  the  decrease  in  effi¬ 
ciency  caused  by  the  DCFC.  When  0<J^Jc  or  J  >  JM,  the  curves 
of  the  power  density  and  efficiency  of  the  hybrid  system  are  over¬ 
lapped  with  that  of  the  sole  DCFC.  When  Jc<J<  Jm,  the  power  den¬ 
sity  and  efficiency  of  the  hybrid  system  are  larger  than  either  that 
of  the  DCFC  or  that  of  the  TEG.  It  clearly  shows  that  the  perfor¬ 
mance  can  be  effectively  enhanced  by  coupling  a  TEG  to  the  DCFC 
for  further  converting  the  waste  heat  into  power  generation  in  spe¬ 
cific  operation  conditions,  and  the  increase  in  the  power  density  is 
more  significantly  than  that  in  the  efficiency.  For  the  typical 
parameters  given  in  Table  1,  the  maximum  power  density  of  the 
hybrid  system  is  318.46  W  m  2,  which  is  about  50.81%  larger  than 
that  of  the  sole  DCFC,  21 1.16  W  m  2.  Meanwhile,  the  efficiencies  of 
the  hybrid  system  and  DCFC  corresponding  to  the  their  maximum 
power  densities  are,  respectively,  38.63%  and  32.30%,  and  the  effi¬ 
ciency  of  the  hybrid  system  is  about  19.60%  larger  than  that  of  the 
sole  DCFC. 

It  is  seen  from  Figs.  4  and  5  that  the  power  density  and  effi¬ 
ciency  of  the  hybrid  system  are  affected  by  the  operating  temper¬ 
ature  of  the  DCFC,  T,  and  the  thermal  conductance  of  the  TEG,  K. 
When  J^Jc  or  J  ^  JM,  the  efficiency  and  power  density  of  the 
hybrid  system  are  increased  with  the  increasing  operating  temper¬ 
ature  T,  while  the  efficiency  and  power  density  are  the  same  as  that 


J  (A/m2) 


Fig.  4.  The  effects  of  operating  temperature  T  on  the  (a)  power  density,  and  (b) 
efficiency  of  the  hybrid  system. 


J  (A/m2) 

Fig.  5.  The  effects  of  heat  conductivity  if  on  the  (a)  power  density,  and  (b)  efficiency 
of  the  hybrid  system,  where  ]s  is  the  stagnation  current  density  that  the  DCFC  . 


of  the  sole  DCFC  for  different  thermal  conductance  K.  The  values  of 
Jc,  Jm  and  JP  shift  to  larger  ones  as  the  operating  temperature  T  or 
thermal  conductance  K  is  increased,  and  JP  is  usually  located 
between  Jc  and  JM.  The  effective  operating  current  density  interval 
A J(=JM  -Jc)  increases  and  moves  backward  as  the  operating  tem¬ 
perature  T  or  the  thermal  conductance  K  is  increased.  It  is  also  seen 
from  Fig.  5(a)  that  there  is  an  optimum  value  of  the  thermal  con¬ 
ductance  K  at  which  the  maximum  power  density  can  be  attained. 
For  the  typical  parameters  given  in  Table  1 ,  the  optimum  value  of  K 
is  found  to  be  between  0.02  W  K  1  m  1  and  0.03  W  K  1  m  Based 
on  numerical  calculations,  more  detailed  values  of  Jc,Jm,  A J,JP  and 
P*Tlilx  under  different  operating  temperature  T  and  thermal  conduc¬ 
tance  K  are  listed  in  Table  2.  Statistical  analysis  may  be  further 
used  to  optimize  the  designing  and  operating  parameters  if  some 
proper  optimization  objectives  and  more  numerical  calculation 
examples  are  given.  Combined  Table  2  and  Figs.  4  and  5,  it  is  inter¬ 
esting  to  note  that  the  power  density  and  efficiency  of  the  hybrid 
system  will  be  still  larger  than  zero  even  when  the  operating  cur¬ 
rent  density  is  larger  than  the  stagnation  current  density  Js  that  the 
DCFC  does  not  contribute  any  more.  For  this  case,  the  DCFC  in  the 
hybrid  system  is  operated  as  a  high-temperature  heat  source  more 
than  an  electrochemical  converter. 

The  influences  of  the  integrated  parameters  c,  and  c2  on  the 
performance  of  the  hybrid  system  are  clearly  shown  in  Fig.  6.  It 
is  observed  from  Fig.  6  that  the  influence  of  c,  is  only  occurred  in 
the  region  of  Jc<J< Jm,  the  power  density  and  efficiency  will 
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Table  2 

The  values  of  jc,  jin,  jp  and  P'mxx  for  different  operating  temperature  T  and  thermal 
conductance  K,  where  ct  =  0.1  W  nr2  K-1  and  c2  =  0.1  W  nr2  K-1. 


T(K) 

K 

(wr1  m- 

Jc 

)  (A/m2) 

Jm 

(A/m2) 

A I-Jm-J 
(A/m2) 

c  Jp 
(A/m2) 

pu 

(W/m2) 

923 

0.01 

530 

888 

358 

762 

292.4 

0.015 

628 

1105 

477 

913 

318.5 

0.02 

719 

1296 

577 

1040 

330.3 

0.025 

803 

1469 

666 

1152 

332.6 

0.03 

882 

1629 

747 

1252 

327.7 

973 

0.01 

582 

974 

392 

832 

322.0 

0.015 

686 

1205 

519 

992 

347.3 

0.02 

782 

1409 

627 

1126 

357.3 

0.025 

870 

1593 

723 

1243 

356.9 

0.03 

954 

1763 

809 

1348 

349.1 

1023 

0.01 

635 

1060 

425 

903 

351.7 

0.015 

745 

1306 

561 

1071 

375.9 

0.02 

845 

1522 

677 

1212 

383.7 

0.025 

938 

1717 

779 

1334 

380.6 

0.03 

1025 

1897 

872 

1443 

369.5 

decrease  as  the  parameter  c,  increases,  while  the  current  densities 
Jc,  Jm ,  and  JP  will  increase  as  the  parameter  c,  is  increased.  When 
the  regenerative  losses  in  the  regenerator  are  negligible,  i.e., 
c,  =  0,  the  fundamental  relationship  between  the  operating  current 


density  of  the  DCFC  and  the  dimensionless  current  of  the  TEG  can 
be  reduced  from  Eq.  (14)  to 

f  - 2ZTj-2Z(T-T0)  +  2f  [n^(l  - iJocrc);- C2(T- To)]  =0.  (20) 

For  this  case,  the  curves  of  the  power  density  and  efficiency 
varying  with  the  operating  current  density  are  shown  by  the  red 
dash  curve  in  Fig.  6. 

Similar  to  the  effect  of  parameter  ci,  the  influence  of  the  inte¬ 
grated  parameter  c2  on  the  performance  of  the  hybrid  system  is 
also  only  occurred  in  the  region  of  Jc<J< Jm,  the  power  density 
and  efficiency  will  decrease  with  the  increasing  parameter  c2, 
while  the  current  densities  Jc,Jm,  and  Jf>  will  increase  as  the  param¬ 
eter  c2  is  increased.  When  the  heat-leak  from  the  DCFC  to  the  envi¬ 
ronment  is  negligible,  i.e.,  c2  =  0,  the  fundamental  relationship 
between  the  operating  current  density  of  the  DCFC  and  the  dimen¬ 
sionless  current  of  the  TEG  can  be  reduced  from  Eq.  (14)  to 

/  - 2ZTj-2Z(T-To)  +  2f  [n^(l  - r,DCFc)J-Ci (T— T0)J  =0.  (21) 

For  this  case,  the  curves  of  the  power  density  and  efficiency 
varying  with  the  operating  current  density  are  shown  by  the  blue 
dot  curve  in  Fig.  6. 

When  both  the  heat-leak  from  the  DCFC  to  the  environment  and 
the  regenerative  losses  in  the  regenerator  are  negligible,  i.e.,  Cj  =  0 
and  c2  =  0,  the  fundamental  relationship  between  the  operating 
current  density  of  the  DCFC  and  the  dimensionless  current  of  the 
TEG  can  be  reduced  from  Eq.  (14)  to 

f  ~  2ZTJ  -  2 Z(T  -r0)  +  2^[n^(l  -  t/DCFC)j]  =  0.  (22) 

For  this  case,  the  curves  of  the  power  density  and  efficiency 
varying  with  the  operating  current  density  are  shown  by  the  black 
solid  curve  in  Fig.  6. 

4.  Conclusions 

A  DCFC-TEG  hybrid  system  mainly  composed  of  a  DCFC,  a 
multi-couple  TEG  and  a  regenerator  has  been  proposed  to  effec¬ 
tively  recover  the  waste  heat  generated  in  the  DCFC.  By  consider¬ 
ing  the  main  thermodynamic-electrochemical  irreversibilities, 
the  analytical  expressions  of  performance  parameters  for  the 
hybrid  system  are  derived,  from  which  the  general  performance 
characteristics  are  revealed.  The  operating  current  density  interval 
of  the  DCFC  that  the  TEG  exerts  its  function  is  determined.  It  is 
found  that  the  performance  of  the  DCFC  can  be  significantly 
improved  by  coupling  a  TEG  to  the  DCFC  for  additional  electricity 
generation,  and  the  increase  in  the  power  density  is  more  signifi¬ 
cantly  than  that  in  the  efficiency.  The  influences  of  the  operating 
current  density  and  operating  temperature  of  the  DCFC,  thermal 
conductance  of  the  TEG,  and  two  integrated  parameters  c,  and  c2 
on  the  performance  of  the  hybrid  system  are  discussed.  The  results 
obtained  in  the  paper  may  offer  some  theoretical  guidance  for  the 
performance  improvement  of  DCFCs  through  cogeneration 
systems. 
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